A Molecular Hydrodynamic Theory of Supercooled Liquids and Colloidal Suspensions 

under Shear 
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We extend the conventional mode-coupling theory of supercooled liquids to systems under station- 
ary shear flow. Starting from generalized fluctuating hydrodynamics, a nonlinear equation for the 
intermediate scattering function is constructed. We evaluate the solution numerically for a model 
of a two dimensional colloidal suspension and find that the structural relaxation time decreases as 
•y~ v with an exponent v < 1, where 7 is the shear rate. The results are in qualitative agreement 
with recent molecular dynamics simulations. We discuss the physical implications of the results. 
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Recently, there has been an explosion of interest in 
understanding out of equilibrium properties in super- 
cooled liquids. In general, the nonequilibrium behav- 
ior of a glassy system is characterized by a violation of 
the fluctuation-dissipation theorem and the absence of 
time translation invariance. In the particular case where 
the system is subjected to a homogeneous, steady shear 
flow, time translation invariance is recovered. This sim- 
pler nonequilibrium situation is of interest for two rea- 
sons. First, understanding the rheological properties of 
complex fluids such as colloidal suspensions and poly- 
mers at a microscopic level is important for the design 
and control of new materials. On a more fundamental 
level, it has recently been suggested that for supercooled 
liquids there are fundamental connections between the 
standard thermodynamic control variables of tempera- 
ture and density in the equilibrium case and steady state 
shear out of equilibrium A major goal of this work is 
to develop a theory that provides an explicit microscopic 
connection between the temperature, density and shear 
rate in a supercooled liquid. 

Dense colloidal suspensions are known to exhibit weak 
shear thinning behavior. Such behavior is predicted for 
simple liquids as well, but the effect is too small to ob- 
serve at temperatures well above the glass transition. For 
supercooled liquids, however, the situation is different. 
Recent numerical simulations have revealed anomalous 
rheological behavior in supercooled liquids. Yamamoto 
and Onuki|2| and Berthicr and Barrat |3j have simu- 
lated supercooled liquids under strong stationary shear 
flow and observed a characteristic shear dependence of 
the structural relaxation time and the shear viscosity, 
T a ,r\ tx A i~ v 1 where r a is the structural relaxation time, 
77 is the shear viscosity, 7 is the shear rate and the ex- 
ponent v is empirically found to range between 2/3 and 
1. An abstract schematic approach based on the exactly 
solvable p-spin spin glass has been proposed and studied 
by Berthier, Barrat and Kurchan 4|. This model predicts 
v = 2/3 in agreement with the lower bound found in the 
simulations of Berthier and Barrat 3J . Since this model is 



schematic, it cannot be used to understand in detail the 
microscopic relationship between fluid structure and dy- 
namics as a function of thermodynamic control variables 
and external driving. 

In this letter, we shall generalize the mode-coupling 
theory developed to describe the fluctuations in an equi- 
librium state to that of a system under a stationary shear 
flow. Our starting point is generalized fluctuating hydro- 
dynamics. Using several approximations, we obtain a 
closed nonlinear equation for the sheared generalization 
of the intermediate scattering function. In this letter, we 
shall consider both normal liquids as well as the over- 
damped Brownian behavior of a colloidal suspension in 
the absence of hydrodynamic interactions. Numerical re- 
sults will only be presented for the Brownian case, but 
the more general results derived here could be used to 
make quantitative contact with recent molecular dynam- 
ics simulations. 

Consider the shear flow given by 



v (r) = r-r= ( 7 y,o,o), 



(1) 



where (T) a p — A /5 ax 5/3 y is the velocity gradient matrix. 
The hydrodynamic fluctuations for density p(r, t) and the 
velocity field v(r, t) obey the following set of Langevin 
equations ["" 
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(2) 



where Co is the collective friction coefficient for colloidal 
particles and ifi(r,t) is the random force. The Co term 
is specific for the colloidal case. In the case of atomic 
liquids, the friction term should be replaced by a stress 
term which is proportional to the gradient of velocity 
field multiplied by the position dependent shear viscosi- 
ties. Both cases, however, lead to the same dynamical 
behavior at long time scales. We neglect the weak r de- 
pendence of the friction coefficient that could arise from 
hydrodynamic interactions between colloidal particles in 
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the Brownian case. The first term in the right hand side 
of the equation for the momentum is the pressure term 
and T is the total free energy in a stationary state. Here 
we assume that the free energy is well approximated by 
that of the equilibrium form and is given by a well-known 
expression; 



I3T ~ jdr p(r){lnp(r) - 1} 

- ^ Jdrt Jdr 2 5p{v 1 )c{r 12 )5p{v 2 ), 



(3) 



where (3 = 1/k^T and c(r) is the direct correlation 
function. Under shear, it is expected that c(r) will 
be distorted and should be replaced by a nonequilib- 
rium, steady state form c non eq(r), which is an anisotropic 
function of r. The effect of shear on the static cor- 
relation functions has been studied^, 0] and is found 
that the distortion is characterized by the Peclet num- 
ber Pe= ja 2 /Dq, where a is the diameter of the particle 
and D$ — k^T/d^o is the diffusion constant. Thus, if the 
Peclet number is small, the above assumption is valid. 
If necessary, the full anisotropic steady state structure 
may be used. By linearizing eq. J2J around the stationary 
state as p = po + 5p and v = v + <5v, where po is the 
average density, we obtain the following equations, 



— - k • r • — ) s Pk (t) = -ikJ k (t), 



* -k.r.JL + k.r-kjA./. 



ik 



m(3S(k) 



Sp^t) 



--L [it- qc(g)5p k _ q (i)5p q (i) - ^J k (i) + f Rk (t), 

TO/3 J q TO 

(4) 

where c(q) is the Fourier transform of c(r), k = k/|k|, 
Jk(t) = pok • £v k (t) is the longitudinal momentum fluc- 
tuation, and J q = J dq/(27r) d for a <i-dimensional sys- 
tem. Note that our approximate equation does not con- 
tain coupling to transverse momentum fluctuations even 
in the presence of shear. 

In order to construct equations for the appropriate 
correlations from the above expressions, an approximate 
symmetry is necessary. In the presence of shear, transla- 
tional invariance is violated. In other words, correlations 
of arbitrary fluctuations, /(r, t) and <?(r, £), do not satisfy 
(/(r,i) 5 (r',0)) ^ (/(r-r',i)t/(O,0)). Instead, we shall 
assume that the following symmetry is valid 0, 



expressed as 

</k(*)flk'(0)> = (/k(t)5 k(t) (Q)} x «5 k(t) , k , 



(/k' ( - t) («0)) xS 



(6) 



k,k'(-t) 



where k(i) = exp[*rt] -k = k + jtk x e y , where *r denotes 
the transpose of T and 5 k . k ' = (2Tr) d V~ 1 5(k — k') for a 
system of volume V. 

Eq. JSJ states that the fluctuations satisfy translational 
invariance in a reference frame flowing with the shear 
contours. This approximation holds for long wavelengths 
where the direct interactions between particles are not 
important. On the other hand, for correlations between 
particles separated by molecular length scales, this is not 
generally true. The validity of this approximation for 
molecular length scales should be systematically exam- 
ined in the future. 

Using this approximation, it is straightforward to con- 
struct the mode-coupling equations for the appropriate 
correlation functions. We shall derive the equation for 
the intermediate scattering function defined by 



F(M) = A(«5 Pk( _ t) (i)^ k (0)), 



(7) 



where N is the total number of the particles in the sys- 
tem. Note that the wave vector in dp^t) is now replaced 
by a time-dependent one k(— t). 
Eq.Q has a nonlinear term 



Rk(t) 



m(3 



i(k • q)c(q)«Sp k _ q (i)^ q (i). (8) 



This term can be renormalized with the definition of a 
generalized friction coefficient following the standardpro- 
cedure of derivation of the mode-coupling equations 0-0]. 
To lowest order in the fluctuations, we have 



dt k 1 dk 



k-r-k j k (t) = — 



mpS(k) 

- f dt' /dk' C(k,k',i-t')./k'(*')+fk(*) ) 



Spk(t) 



(9) 



where £(k, k', t) is the generalized friction coefficient and 
ffl k (£) is a corresponding random force. £(k, k', t) is given 
by the sum of the bare friction coefficient and the mode- 
coupling term as 

C(k, k', t) = Co x 26(t) + S((k, k', t) (10) 

with the mode-coupling contribution given by 



*C(k,k',t) = 



m 2 j3 
N 



(R k (t)Rl(0)). 



(11) 



(/(r, t)g(r', 0)) = (/(r - r'(t),t)g(0, 0)), (5) Substituting eq.© into the above expression, we obtain 



where we defined the time-dependent position vector by 
r(t) = exp[rt] • r = r + jtye x , where e x is an unit vector 
oriented along the cc-axis. In wave vector space, this is 



S((K k'i) = / / k • qc(q)k' • q'c(q') 

PN J q J q , (12) 

x {5p^{t)5 Pci {t)5pl,_AQ)5p*,(Q)). 
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This involves a four point correlation function. Using the 
Gaussian approximation, this can be decomposed into a 
product of two-point correlation functions as 

(<Jpk-q(*)W*)*Pk'-q' WW )) 

~iV 2 F(k(t)-q(t),t)F(q(t),t) (13) 

x {<5k(t),k'^q(t),q' + <5ik(t),k'£q(t),k'-q'} j 

where use has been made of the translational invariance, 
eq.©. 

Substituting ca. l|13() back to ea. l|12ll . we obtain 



S((k,k',t) = S((k,t)S W)x 



(14) 



with 



6((k, t) = ^ y {k • qc(q) + k • (k - q) c(k - q)} 

x k(t) • q(t)c(q(t))F(k(t) - q(t), i)F(q(t), t) 
= g| V(k,q)V(k(t),q(t))F(k(t)-q(t),t)F(q(t),t), 

(15) 

where V(k, q) is the vertex function given by 

V(k, q) = k • {qc(q) + (k - q) c(k - q)} . (16) 

From these results and eq.|0, the equation for the cor- 
relation function, 



C(M) = -<J k( _ t) (t)n£(0)) 



is given by 
dC(k, t) 



dt 



k(-t) • r • k(-t)C(k, t) — Stt^-tt^^ *) 



m/3S(k{-t)) 
1 /"* 

— / dt' $C(k(-t),t-t )C(k,t'). 
to Jo 



(17) 



Note that in the above equation, the differential operator 
k ■ r • d/dk disappears because 

dc(k,*)_«7(k,t)_ >r a (lg) 



dt 



at 



flk(-t) 



Likewise, the continuity equation (the first term in 
eq.Q) can be written as 



dF(k, t) 
dt 



-ifc(-t)C(k,f). 



(19) 



This equation together with eq. (|1 T|) comprises the closed 
set of the mode-coupling equations for F(k, t) and C(k, t) 
under shear. 



For colloidal suspensions the relaxation time of the mo- 
mentum fluctuations is of the order of r m = to/Co and is 
much shorter than the relaxation time for density fluctu- 
ations which is of the order of or longer than Td = a 2 /Do. 
For the time scale of interest, k(— t) • T • k(— t) as well as 
the inertial term can be neglected in the equation for the 
momentum fluctuations since jT m <C 1 in realistic situa- 
tions. Thus, the equation for the momentum fluctuations 
may be written as 



- f dt' (5C(k(-t),t-t')C(k,t'). 
Jo 

Substituting this back into ea. H19(l . we arrive at 
dF(k,t) D k(-t) 2 



(20) 



dt 



S(k(-t)) 



F(k,t) 



- ( dt' M(k(-t),t-t') 
Jo 



, .dF(k,f) 



(21) 



dt' 



where 



M(k,t)=^0-A y / q V(k ) q)V(k(t),q(t)) 

xF(k(t)-q(t),t)F(q(t),t). 

Eqs.J2U an( i are the major result of this letter. In 
the absence of the shear, they reduce to the conventional 
mode-coupling equations 

In order to study the shear thinning effect in the su- 
percooled state, we shall solve eas. H21|) and (1221) numeri- 
cally. Solving this equation is more difficult than solving 
the corresponding equation in the equilibrium case be- 
cause the wave vectors are distorted by shear and the 
system is not isotropic. For simplicity, we shall con- 
sider a hypothetical two-dimensional colloidal suspension 
which is simple to handle numerically but still undergoes 
an ergodic-nonergodic transition below a certain density. 
The shear flow occurs in the x direction. We have chosen 
the following form of the static structure factor S(k): 



S(k) = S PY (k + k 0l apa)f(k - k c ), 



(23) 



where SpY(k,p) is the static structure factor for a hard- 
sphere system at the density p obtained from the Percus- 
Yevick closure, ko and a are parameters which were cho- 
sen in such a way that S(k) is short-ranged and has 
broader peak. f(k — k c ) is a cut-off function which 
makes S(k) approach unity smoothly for wave vectors 
larger than the cut-off k c . The choice of S(k) mimics the 
shape of S(k) of real systems although it does not sat- 
isfy sum-rule restrictions. In our calculation, we chose 
ko = 4.0, a = 32 and k c = 4.0. For this system, the 
ergodic-nonergodic transition occurs around a "density" 
p c a 2 ~ 1.2 x 10~ 2 in the absence of shear. In Figure 1, 
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we show the behavior of F(k,t) for pa 2 — 1.15 x 10~ 2 , 
slightly below p c . The wave vector is k = (0,2a). Since 
k x = 0, the expression for F(k, t) is equivalent to that 
in the quiescent state, F oq (k, t) = N~ 1 {6pk(t)6p^(0)). 
Thus, there is no direct effect from the convection term 
but, due to the nonlinear coupling through the mode- 
coupling term M(k, t), a strong shear dependence of the 
relaxation time can be seen. The dependence of the 
structural relaxation time T a ("f) on the shear rate 7 is 
estimated from the value where F(k, r a ) = e _1 . For the 
particular case of pa 2 = 1.15 x 10 -2 we find the power 
law r a ~ 7 — " with v ~ 0.8 for Pe> 10 -3 and r Q sat- 
urates to the equilibrium value at Pe< 10~ 3 . This is 
similar to the behavior reported in recent molecular dy- 
namics simulations. For values of pa 2 > 1.2 x 10~ 2 (not 
shown) , we find that the exponent v saturates at a higher 
value, in agreement with the simulations of Bcrthier and 
Barrat 0. 

The physical picture that emerges from the molecu- 
lar hydrodynamic theory developed here is simple. The 
shear flow perturbs and randomizes the coupling between 
different modes. Physically, this perturbation dissipates 
the cage that transiently immobilizes particles. Mathe- 
matically, this is reflected through the time dependence 
of the vertex, which vanishes as t — > 00. This simple 
picture illustrates the essence of the mode-coupling ap- 
proach to the shear thinning effect in simple supercooled 
systems. Note that even for fluctuations orthogonal to 
the direction of flow, thinning occurs due to the coupling 
of fluctuations in all directions. In this sense, the picture 
of cage breakup in a supercooled liquid due to external 
flow is quite different from that of dynamic critical phe- 
nomena under shear, in which the faster relaxation occurs 
solely because the fluctuations are stretched out by the 
shear flow and pushed to larger wave vectors where faster 
relaxation occurs. 

In this letter, we have derived an approximate mode- 
coupling theory for a supercooled liquid under steady 
shear flow. The most important assumption is the use 
of approximate translational invariance, eq.JHJl- This al- 
lows one to derive a nonlinear integro-differential equa- 
tion for F"(k, i) similar to the one for the equilibrium 
state. The numerical analysis for a hypothetical two- 
dimensional colloidal suspension has been carried out and 
a typical behavior of F(k, t) was shown to be consistent 
with recent simulations. The relaxation time is found 
to have the strong shear dependence. More systematic 
and thorough analysis of the numerical solutions of the 
mode-coupling equations are left for future work^lj. 
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FIGURE CAPTIONS 



FIG. 1: Normalized F(k,t) for k = (0,2a) for various shear 
rates 7. The "density" is pa 2 = 1.15 x 10~ 2 . From the right 
to the left, Pe = ia 2 /D = 0, 1CT 4 , 1CT 3 , 1CT 2 , HT 1 , and 1. 
The results for Pe= and 10 -4 are almost indistinguishable. 
The time t is scaled by a 2 /Do. 



